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Abstract — This paper extends the consensus framework, widely 
studied in the literature on distributed computing and control 
algorithms, to networks of quantum systems. We define consensus 
situations on the basis of invariance and symmetry properties, 
finding four different probabilistic generalizations of classical 
consensus states. We then extend the gossip consensus algorithm 
to the quantum setting and prove its convergence properties, 
showing how it converges to symmetric states while preserving 
the expectation of permutation-invariant global observables. 

Index Terms — Consensus, Quantum Information and Control, 
Markov Processes, Distributed Algorithms. 



I. Introduction 

Among the recent trends in control and systems theory, the 
field of distributed control, estimation and optimization on 
networks has stimulated an impressive amount of research, 
see e.g. 0, Q, 0, 0, 0, 0, 0, 0, @. A basic task 
for distributed information processing is reaching consensus 
on some common control objective, shared value or slack 
variable. The present paper extends this well-studied consen- 
sus problem (see e.g. (TOj, flj], H3, [13]) to networks of 
quantum systems, a special case of which would be classical 
probability distributions. 

Exploring the links between information processing tasks 
and stochastic dynamics on networks has recently opened new 
research directions towards "distributed" quantum information 
applications. Among these, we recall quantum computation 
03), [15 1 in its potential implementation via dissipative means 
[16 1, and its connection to quantum random walks [17|, |18|. 
Other specific tasks include noise protection and dynamical 
error-correction 021, EH, ED, E2, E3, entanglement gen- 
eration through stabilizing dissipative dynamics [24], ||25ll . the 
creation of quantum "gadgets" driven by dissipation l26l . as 
well as most tasks in the control of open quantum systems 

E2, EH, E21, ED, ED. 

In this spirit, we here develop a framework for addressing 
quantum consensus problems: we identify and characterize a 
hierarchy of quantum consensus situations and study how these 
can be reached by suitable dissipative quantum dynamics, 
while preserving some global information on the network. Our 
work offers not only a formal generalization of the well-known 
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classical consensus problem, along with some intriguing in- 
sights on the latter, but also a potentially new viewpoint on a 
number of issues in quantum information. More specifically, 
it ties the structure of symmetric states and correlations [32|, 
ll33ll to their potential generation via quasi-local resources [ 25 1, 
and to the relaxation of many-body systems to thermal states. 
While this work is mainly meant to provide solid theoretical 
foundation to quantum consensus problems, links between the 
proposed theory and potential applications will be outlined in 
the conclusions. As an interesting by-product, the present work 
also provides a viewpoint on consensus for classical systems in 
the context of probability distributions instead of deterministic 
states. 

An attempt to lift the consensus problem to the quantum 
domain has been presented in J34). It is based on a "cone 
geometry" approach, viewing quantum Kraus maps as the non- 
commutative generalization of Markov chains that model con- 
sensus algorithms. The authors show how Birkhoff's Theorem 
and Hilbert's projective metric lead to a general convergence 
result and contraction ratio estimation. However, by describing 
the dynamics of the whole system of interest as governed by a 
single Markov transition mechanism, this formulation does not 
account for subsystem structure or network connections in the 
quantum setting. It therefore defines consensus as asymptotic 
convergence to a scalar multiple of the identity: for quantum 
states this corresponds to a fully mixed, most uncertain state 
which is not the typical desired target for quantum information 
processing applications. 

In the present paper, we approach quantum consensus from 
an "operational", multi-agent control perspective, starting from 
the basic classical ingredients: a network of subsystems, an 
interaction protocol with locality constraints, and a target 
consensus situation. As a first step, Section [II] provides four 
possible ways to generalize the concept of a consensus state to 
the quantum domain and explores their connections, establish- 
ing a hierarchy of quantum consensus definitions. Section III 
presents the quantum open-system dynamics and the locality 
notions that we employ to describe the interactions between 
the quantum "agents". Typical methods for interaction selec- 
tion and timing are also introduced, in analogy with classical 
consensus. This provides all the ingredients needed to specify 
the general properties of an evolution that achieves quantum 
average consensus. The symmetry-based reformulations at 
each step provide an alternative interpretation of classical 



consensus, and in Section IV this leads us to a quantum 
generalization of the classical gossip algorithm. We prove its 
convergence to symmetric-state consensus while preserving 
the expectation of any permutation-invariant observable. We 
further show how the algorithm can be explicitly seen as a 
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generalization of the classical one and the classical conver- 
gence results can be used to prove a weaker convergence 
property. The section concludes with an example of the gossip 
algorithm working on a four-qubit network. In Section |V| we 
summarize the main results and provide an outlook on possible 
developments and applications for our framework. A small 
tutorial on quantum systems modeling is given in Appendix 
[A] along with an overview of the notations and conventions 
we use. 

II. Consensus States 

A classical consensus state for a multipartite system is 
one in which the states of all the subsystems, often called 
agents, are the same (although not necessarily stationary). 
Since in general a quantum state for a multipartite system 
cannot be factorized into subsystem states (see Appendix [A), 
the definition of consensus must be carefully reconsidered. We 
here present an operational approach that leads to a hierarchy 
of definitions for what can be claimed to be "quantum con- 
sensus". 

A. Defining Classical Consensus 

The consensus problem for classical systems is most 
typically formulated along the following lines iflOl . |12|. 
Consider m subsystems, each one associated to a state 
given by a configuration variable Xk E R". These sub- 
systems evolve through bilateral interactions, according to 



some networking scheme and dynamics (see Section IV i, and 
reach a consensus state if they converge to the set C = 
{ (xi, X2, x m ) : Xj = XfcVj, k}. Furthermore, the agents 
are said to compute an average consensus if for given initial 
states (xi(0), ^2(0), x m (0)) they converge to the particular 
equilibrium (x, x, x) E C where x = A x i(fy- 

Alternatively, consensus can be characterized as invariance 
w.r.t. subsystem permutations. Let *}3 denote the set of all 
subsystem permutation operators, i.e. each P^ E *p is asso- 
ciated to some permutation tt of the integers 1,2, ...m such 
that P 7T (x 1 ,x 2 ,...,x m ) = (a; 7r (i),a; 7r (2),... ) ^( m )) for any 
xi,X2, ■■■,x m . Denoting the joint state of the Xk E M™ by a 
vector x € M. mn , each P n E can be written as an ran x mn 
matrix resulting from the Kronecker product of some m x m 
permutation matrix with the n x n identity matrix. Then we 
can define consensus as: 



C = {xE 



P w x = x for all Pit E ^P} . 



(1) 



It is standard that checking P^ x — x for all pairwise 
permutations P v is sufficient to guarantee P^x = x for all 

Average consensus is a fundamental element of distributed 
computation because for each linear function Q : M. mn — > K s 
which is invariant under all subsystem permutations, there 
exists q : W l — > R s such that Qx — qx. Thus such 
Q can be evaluated locally by each agent once they have 
reached average consensus. The fact that all permutations can 
be obtained by concatenating pairwise permutations, already 
suggests that average consensus is computable with standard 
distributed algorithms. 



B. Quantum Consensus Definitions and their Relationships 

Defining what a consensus situation ought to be in a 
quantum "network" is not a straightforward task. More than 
one definition may be appropriate depending on the type of 
symmetry we are seeking. Following the analogy with the 
classical case can help, but quantum measurement outcomes 
are intrinsically stochastic, so we must consider probabilistic 
consensus situations from the beginning. Let us explore 
different options by first discussing a simple case. 

Example 1: When is a quantum network in consensus? 

Consider a multipartite quantum system composed of three 
qubits, with associated Hilbert space T~L = C 2 <8> C 2 <g> C 2 , and 
three observables of the form = a z ®I®I, — I®a z ® 
I, cj( 3 ) = I®I®cr Zl where the Pauli matrix a z = diag(l, — 
These correspond to observables of the quantity associated to 
a z for each of the subsystems. It seems natural to say that the 
system is in consensus with respect to the expectation of a z 
if 

Tr(pa (1) ) = Tr(pa^) = Tr(pcr^). (2) 

The conditions for this to happen can be worked out explicitly 
in terms of the diagonal elements of the state p. In particular it 
is easy to check that all the following states satisfy condition 
(0: " 

/ = i/O(|0) + |l))«0| + <l|)<»(|0) + |l))«0| + <l|); 
P B = ~/®(|0,0> + |l,l))(<0,0| + (l,l|); 



(|0,0,0)<0,0,0| + |1,1,1)(1,1,1|); 



p E = 10,0,0X0,0,01; 

P F = ^(|0,0,0} + |1,1,1»«0,0,0| + (1,1,1|). 

All these states, except p E , have Tr(pa^) — for i = 1, 2, 3. 

The above requirement can be strengthened by requesting 
(|2| to hold when a z is replaced by any observable a E *B(C 2 ) 
in the definition of This is equivalent to 

imposing that the reduced states for the three subsystems are 
the same (a formal proof is provided in the next subsection). 
It is then easy to check that p B , p c , p D , p E , p F satisfy this 
requirement, while p A does not. In fact the reduced states for 
p A are: 

2 7 < ^ = "3='( 

In the light of ([TJ, another potential definition of quantum 
consensus would require the state to be symmetric, i.e. invari- 
ant under any permutations of the subsystems. This choice can 
be motivated by the classical case, where the consensus state 

'A qubit is a quantum system associated to a two-dimensional Hilbert 
space W ~ C 2 ; a standard basis for the latter is conventionally given by the 
vectors {|0), |1)} C H. The traceless unitary and hermitian Pauli operators 
<T x ,a y ,a z , completed with the identity operator, form an orthonormal basis 
for the operators on H. Explicitly, a x = |1)(0| + lOXll.o-y = i|l>(0| - 
i[0)(l|,<T. = [0)<0[-|l)(l[. 



p? = ~i, P2=pi = aio) + |i))«o| + (i|). 
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is indeed permutation invariant. Among the states defined in 
Example 1, only p c ,p D ,p E ,p F are permutation invariant. 

Lastly, one might desire to have subsystem agreement 
not only on the observable averages, but also on each 
measurement. Namely, we want that each projective 
measurement of the (commuting and hence compatible) 
observables a\ , 02, 03 gives perfectly correlated results 
for the three subsystems. That is, among all possible 
measurement results {— l,+l} x3 , only (—1,-1,-1) and 
(+1,+1, +1) have a nonzero probability to occuiH For 
instance if the system is in an entangled state of the 
form p = HO, 0,0) + j3\l, 1, l))(a*<0, 0, 0| + /3* (1, 1, 1|), 
with \a\ 2 + \f3\ 2 = 1, the outcome of a measurement 
operation is (—1,-1,-1) with probability \a\ 2 and 
(+1,+1,+1) with probability |/3| 2 . With a = /3 = i 
we get p F above. The states p A , p B and p c do not 
satisfy this definition of consensus; indeed, for these three 
states, the distribution of measurement results for qubit 
1 is independent of the measurements on the other two 
qubits, e.g. (+1,-1,-1) and (—1,-1,-1) have the same 
nonzero probability so uncorrected measurement results 
can happen. Note however that also mixed states can 
lead to correlated results, if they express perfect classical 
correlations, as is the case for any state of the form 
p G = p|0,0...0)<0,0...0| + (1 -p)|l,l...l)(l,l...l| 
with < p < 1. For three qubits with p — h we get p D 
above. □ 

Let us formalize the ideas emerging from the former 
example. Consider a multipartite system composed of m 
isomorphic subsystems, labeled with indices i = l,...,m, 
with associated Hilbert space H m := Hi ® • • ■ ®U m ~ U® m , 
with dim('H J :) = dim('H) = n and n ^ 2. We shall refer 
to this multipartite system as to our quantum network. For 
any operator X e S("H), we will denote by X® m the tensor 
product X (g) X (£)...(£> X with m factors. Given an operator 
a g ^{Ti), we denote by a^ 1 ' the local operator: 



r (i) := 



j®(m— i) 



Permutations of quantum subsystems are expressed by a 
unitary operator € il(H), which is uniquely defined by 



u v {x l ®...®x m )ul = x <1) 



X 



7r(m) 



for any operators X\,... X m in 03 ("H), where 7r is a permuta- 
tion of the first m integers. A state or observable is said to be 
permutation invariant if it commutes with all the subsystem 
permutations. It is worth noting that given any observable 
Q £ $)(H m ), we can define a permutation invariant observable 
X by considering: 



X = 



1 



UlQU„ 



(3) 



2 The set {ci, C2, C3, ...} xn is the cartesian product of {cj, 02,03, ...} 
by itself n times, i.e. the set of n-tuples with components taken from 
{ci,c 2 ,c 3 , ...}. 



Definition 1 (aEC): Given <r e $ 8('H), a state p e D{H m ) 
is in a -Expectation Consensus (crEC) if: 

Tr(o-Wp) = . . . = Tr(a (fe) p)- 

Definition 2 (RSC): A state p g D("H m ) is in Reduced 
State Consensus (RSC) if, defining the reduced states = 

Tr (®^) (p) ' We haVe 

Pi = • • • = Pm- 

Definition 3 (SSC): A state p € T)(H m ) is in Symmetric 
State Consensus (SSC) if for each unitary permutation we 
have 

u, P ul = p 

Definition 4 (aSMC): Given an observable a with spectral 
decomposition a = SjHj € T)(7i), all different, a state 
/? G 35 ("H" 1 ) is in Single a -Measurement Consensus (crSMC) 
if: 



Tr(nf nf = Tr(nf V), 



(4) 



for all k, I £ {1, . . . , m}, and for each j. 
The definition of crSMC requires that the outcomes of a 
measurements on different subsystems be exactly the same for 
each trial. Indeed, in this last definition, the right-hand side of 
Q is the probability of obtaining Sj as a measurement result 
on both subsystem £ and k (note that IX^ and 11^ commute, 
so this joint measurement nj fc ^II^ is well-defined). If those 
two probabilities are equal, then necessarily the probability of 
Sj on k conditional to observing sj on I is one (assuming that 
Jl^pllf ^ 0; that special case is trivial and can be treated 
separately). 

One may wonder why, as we did for the consensus in 
expectation, we did not try to strengthen the crSMC property, 
requiring perfect correlations of outcomes for any local mea- 
surements. In fact it can be shown that states satisfying this 
property do not exist. 

Proposition 1: For any finite-dimensional quantum system, 
there does not exist a state p that satisfies crSMC for all a. 
The proof makes use of Proposition |4]below, thus we postpone 
it to Appendix [C] □ 

Remark: It is worth remarking how all these definitions 
could be given for classical systems, in the context of consen- 
sus for random variables or for probability distributions of the 
state values. In this case, for example, aEC would require the 
expectation of a set of random variables, each one associated 
to a subsystem, to be the same in all subsystems; RSC would 
require the marginal distributions on each subsystem to be 
equal; and SSC would require that the joint probability dis- 
tribution is invariant with respect to subsystem permutations. 
In some sense the definition of ctSMC is the closest to the 
classical case, as it requires perfect agreement on the outcome 
of a set of random variables for each measurement. 

All the states in our example satisfy tr z EC, all but p A satisfy 
RSC, p c to p F satisfy SSC, and p D to p F satisfy cr z SMC. 
There obviously seems to be a hierarchy in these definitions, 
and the following properties are meant to better characterize 
them. We shall use the notation Cxyz C i D('H) to indicate the 
set of states that satisfy the properties of the acronym XYZ, 
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e.g. Crsc is the set of a ll Reduced State Consensus states. Note 
that, since all consensus definitions involve linear constraints, 
the consensus sets are all convex. 

Theorem 1: The following chain of implications holds: 



SSC 



RSC 



erEC, 



while the converse implications are not true in general. 

Proof: SSC =>■ RSC: If U^pUl = p for each permuta- 
tion, consider in particular Uf^f.) that swaps subsystems I and 
k. Then 

Pk = Tr 0j ^ ^ (p) = Tr 0j ^ K . {U {m pU ] {m ) = p t , 

and the reasoning can be repeated for any pair. RSC => OSC 
is immediate by definition. States p B and p A from Example 
1 provide counterexamples for the converse of the first and of 
the second implication, respectively. □ 

The converse result for the relation between RSC and erEC 
goes as follows. 

Proposition 2: A state is RSC if and only if it is erEC for 
all ae 33(H). 

Proof: RSC implies crEC for all cr by definition. Con- 
versely, assume p to be erEC for any a. The reduced state for 
any subsystem can be written as pk — 2~2e TMpcr^)^, where 
{ae} is an orthogonal, hermitian basis of operators for 05(H). 
Since by hypothesis Ti(pa^) does not depend on k, then all 
the reduced states are identical. □ 

A converse result for the relation between RSC and SSC 
cannot be set up in general. Instead, we have the following 
property. 

Proposition 3: RSC states with pk a pure state for each k 
are also SSC states. 

Proof: If pk is a pure state for each k, then necessarily 
p = \ip){ip\ with \ip) = \ipi) ® \tp 2 ) ® •■• ® \ip m ) f° r some 
|^fc), k = 1,2,..., to. If in addition we require RSC, then 
we need \i/)k){i>k\ = Pk = Pj = \ipj){ipj\ for all j,k, thus 
|"0fe) = IV'j) U P to an irrelevant phase factor for all j, k and 
particle permutation indeed leaves \tp) invariant. □ 

The structure of the crSMC states is further characterized in 
the following proposition. 

Proposition 4: A state is in crSMC if and only if, defining 



n. 



■sym 



^1 7 



it holds 



Tr(II sym( o) = 1, 



or equivalently 



n sym pn 



sym 



flsymP P' 



(5) 



(6) 



Proof: Note that the properties Tr(n sym/ o) = 1 and 



n s ym/°n sym = II sym p = p are equivalent because II sym is 
an orthonormal projector and p is self-adjoint positive semi- 
definite with unit trace. Assume |5) to hold. Along with the 



identities IL? ] U sym = nf ] Uf m = IT 



= Tr(n (£ 

= Tr(n (£ 

= Tr(ITy 

= Tr(nf 



this gives: 



nf } n symP ) 
nf m P ) 



n s ym/^) 

p), 



for all j, k,£. Hence, the SMC definition Q indeed holds. 
On the other hand, suppose that |5]l does not hold. This means 
that Tr((7 — IT sym )p) > 0. We want to show that this implies 



Tr(ITf } nf p) ^ Tr(ny V) 

for some j, k,L Let us write 

i n s vm — 



31, 



E 

.,j m except = 



n, 



n, 



Since Tr((J - II sym )p) > implies that T^IL^ ® ... ® 
Hj m p) > for at least one of the terms in the above sum, 
let us take one such term, denote the corresponding indices as 
{j s } and denote by k,l two subsystems such that jk ^ ji in 
that term. Now writing IL^ £>,... 
where all factors commute, we have: 



®n 5 = n«n( 2) ..'.n> ) 



32 



Tr(m J m > P ) > Tr(m ±; n 



(1) n (2) 

3i h 



n {m) P )>o. 

3m r ' 



Exploiting basic projection properties we thus get: 



Tr(Ili fe) I#p) 



< Tr((l-lf>)irfp) 



31 



Tr(n^p) - Tr(n^^n^V) 

v 3t r/ v 3k 31 r/ 



< Tr(n^p) -Tr(iF (1) iF (2) ...n 



(m) 
3 



p) 



□ 



To conclude, the following proposition connects crSMC to 
the other properties. 

Proposition 5: (a) crSMC implies erEC, while the converse 
is not true. 

(b) ctSMC for a a with nondegenerate spectrum implies RSC, 
while the converse is not true. 

(c) crSMC for a a with nondegenerate spectrum implies SSC, 
while the converse is not true. 

Proof: For (a), we have since Q holds for all k,£: 



Tr(l4 fc) p) = Tr(nf ] uf ] p) 



3 < > \ 3 

By linearity, we thus have: 

Tr(crWp) 



Tr(n^p) 



3 



Sj -Tr(nfp)=£ S ,Tr(Ilfp) 



= Tv(a^p). 

A counterexample for the converse is state p A in Example 1. 

Counterexamples for the converse of (b) and (c) are re- 
spectively states p B and p c in Example 1. For the direct 
statements, given Proposition [T] we know that if (c) is true, 
then (b) must be true as well. Let us then focus on (c). Take the 
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representation of p in the basis associated to a = Y^j s j\j){j\> 
where thus \j)(j\ = IF,, that reads 

P= J! r ji,j2,--^,\ji,j2,--;j m )(ki,k 2 ,...,k m \. 

k ^ , k% , - ■ ■ k rn 

31 > 32> •••3m > 
k ^ , k^ ■ ■ ■ - kjyi 

From Proposition [4j the condition for erSMC writes 

k,j 

so |7]) must reduce to 

P = ^2 Pkj\kk...k){jj (7) 

for some p k j E C. It is straightforward to see that a p of this 
form satisfies SSC, since any element in the sum is invariant 
w.r.t. subsystem permutations. □ 

C. On detecting consensus 

In the quantum setting, the relation between state and 
accessible information is stochastic at a fundamental level. 
In fact, even a "deterministic", maximal information state 
yields probabilistic outcomes for certain observables. As we 
are compelled to use probabilistic notions, consensus can only 
be inferred from stochastic measurement records, and different 
types of consensus require different types of measurement 
statistics. 

The ctEC, requiring only equal expectations for a particular 
observable a on the different subsystems, is trivial to relate 
to measurement results; in particular, checking uEC does not 
require to monitor correlations between measurement results 
on different subsystems. Checking RSC requires statistics for 
a basis of observables for each subsystem; as for crEC, po- 
tential correlations between measurement results on different 
subsystems play no role. On the other hand, distinguishing 
SSC from RSC does require to inspect correlations between 
measurement outcomes at different subsystems - except for 
the case mentioned in Proposition [3] where SSC can be 
inferred directly from the special form of the reduced state 
of each subsystem. A proof of this is provided in Appendix 
[5] For instance, considering the state p B of Example 1, 
measurements of a z on the three subsystems would quickly 
show that the results on subsystems 2 and 3 are always 
perfectly correlated, and show no correlation at all with the 
results on the first subsystem. This difference in correlations 
rules out p B as a candidate for SSC. The definition of ctSMC 
is all about correlations between measurement outcomes at 
different subsystems: the latter must be fully correlated for a 
particular observable a. 

Positively detecting states in Cssc\Co-smCi however, appears 
to be less obvious (except through full state tomography). 
Indeed, a priori it would require to check how permutations 
affect the full state, including not only classical correlations 
but also entanglement (see Appendix [A] for a quick review), an 
intrinsically quantum type of correlation that plays a central 
role in quantum information theory and its applications |14|, 

ma. 



III. Quantum Evolutions for Consensus 

So far we have been concerned with discussing "static" 
properties of consensus states. However, the core of the prob- 
lem is the design of dynamical interactions (or algorithms) that 
drive the system to consensus. In the following sections, we 
will establish which dynamics are needed to drive an arbitrary 
initial state towards a consensus state. In addition to this, as 
in the classical consensus problems, we shall require the final 
state to preserve or express some property dependent on the 
initial state. For example, when classical average consensus 
is reached, each agent locally contains ("has computed") the 
average of the initial state values, which is a global property. 
In the quantum case as well, a goal when reaching consensus 
would often be to retrieve, in the final state of any local 
subsystem, some global information about the initial state. 

A. Classical Dynamics and Locality 

For classical consensus, the starting point is a first-order 
integrator dynamics for each individual agent, of the type: 

x k (t + 1) = x k (t) + u k (t) or f t x k {t) = u k (t), (8) 

for k = 1, 2, m. In this paper we focus on the discrete-time 
case. The inputs u k {t) for the agents can take different forms, 
but they are all based on local information. A notion of locality 
is traditionally introduced by specifying a graph G(V,E) 
whose vertices are the agents 1,2, ...,m and where an edge 
(j, k) E E (unordered pair of vertices) indicates that agents j 
and k are neighbors. Then for each k, input u k is restricted to 
depend only on x k and on the states Xj of agents j for which 
(j, k) E E. Mostly, u k is written as the sum of contributions 
from different edges, u k {t) = J2j-.(j,k)£E fj,k{ x k{t),Xj{t)). 
Directed graphs can be used to model directed information 
flow, i.e. E contains directed edges (ordered pairs of vertices) 
such that (j, k) E E indicates that u k depends on Xj, but not 
necessarily uj on x k (the latter being governed by another 
potential edge, (k,j) E E). Weights Wu jk \ can be associated 
to the edges to model their different relative strengths. Besides 
the given asymmetry through edge weights, one usually forbids 
u k to explicitly use agent identifiers and apply a different treat- 
ment to information coming from or going towards different 
neighbors. This implies 

u k{t) = T,j:(j,k)eE w U,k) f(Xj,X k ) , (9) 

where now / must be independent of j, k. Finally, one may 
request that the influence of agent j on agent k be exactly 
equivalent to the reverse influence, of agent k on agent j. 
Then u k of the type (|9]l features symmetric weights, ^ — 
w (k,j) — leading to an undirected weighted graph — and the 
antisymmetry property f(x, y) — —f(y, x) for any x,y E M. n . 
Such evolution would preserve the average of the state values. 

More generally, locality can be defined on the basis of 
quasi-local operators. Instead of considering a graph, we 
define a set of neighborhoods A/} C {1, ...,m} for j = 
1,...,M, and a quasi-local operator is one that leaves the 
states of all subsystems unchanged except those of one neigh- 
borhood JVj. Then a simple dynamics with local coupling 
would write 

x(t + l)= EjiiV*(*(*)) do) 
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where V; 



for each j is a quasi-local operator 



acting on the neighborhood Afj. Treating all involved agents 
equivalently can be formulated as requiring that for each j, 
the quasi-local operator satisfies PVj — VjP for all the 
permutations P that only permute agents of Afj. The classical 
consensus algorithm ([8|,(|9]l is obtained by taking two-agent 
neighborhoods only, and identifying each Afj with an edge of 
the graph. 



B. Quantum Dynamics and Locality 

According to Schrodinger's equation, isolated quantum sys- 
tems evolve unitarily f36l . 03). However, unitary dynam- 
ics are not enough when we are interested in studying or 
engineering convergence features for a quantum system. A 
more general framework that includes (Markovian) open- 
system evolutions is offered by quantum channels ll37l . 1141 . 
that is, linear, completely positive (CP) and trace preserv- 
ing (TP) maps from density operators to density operators 
£ : T){H m ) -t D(W m ). It can be shown that such maps 
admit an operator sum representation (OSR), also known as 
Kraus decomposition: 



K 



£{ P ) = ]T A kP A 



with 



^4^ 



fc=i 



k=l 



I 



(ID 



where K ^ (dim('H)) 2 . The representation is not unique, but 
all the representations can be obtained as unitary combination 
of the operators of a given one (see [14, Theorem 8.2]). A 
CPTP map is said unital if £ (I) = /. Unital quantum channels 
can always be represented as given by random unitaries [38 1. 
A channel belongs to this class when it admits an OSR with 
operators A k — y/p~kUk, with Uu € il(% m ) and pt > such 
that J2kPk = 1: 

£{ P ) = Y,P^P u l 

k 

Such a map can be thought of as a probabilistic mixture of 
unitary evolutions. 

Given a CPTP map £, we can define its dual map with 
respect to the Hilbert-Schmidt inner product £* : 93(%) — > 
93 ("H) through the relation: 



Tv[A£(p)]=Tr[£HA)p] 



(12) 



This dual map is still linear and completely positive, while the 
fact that £ is trace preserving implies that £t is always unital. 
Considering the dynamics in the dual picture, i.e. with time- 
invariant states and maps acting on the observables, is called 
Heisenberg's picture in the physics literature and provides an 
equivalent description of quantum system evolution. 

We now introduce locality notions for the quantum network. 



Consider the multipartite system introduced in Section II-B 
following l25l . we say that an operator in *B(H) is quasi- 
local if it acts non-trivially only on one neighborhood Afj C 
{l,...,m}: 

Definition 5 (Quantum quasi-local operator): An operator 
V is quasi-local with respect to a set of neighborhoods 



{Afj ,j = 1,2,..., M}, 
{1, 2, M} such that: 

V = 



if and only if there exists j £ 



(13) 



where, with a slight abuse of notation, Vj\f. accounts for the 



non trivial action on Hm and I 



®k 



k^Mj 



C. Timing of operations and evolution types 

In classical consensus, an important aspect is that the graph 
(and the related interaction law) can be time-varying. For 
instance one can assume that all edges are activated for 
the whole time {synchronous update), at the other extreme 
that they are activated one at a time, or some at each time 
{asynchronous update), according to some predefined time- 
varying sequence or by random selection of edges. Interest- 
ingly, convergence properties for all these cases can be linked 
to the connectedness of the "average graph" fOl . 

In the quantum case also this distinction can be made. The 
elementary dynamical interaction that we consider, replacing 
"one edge" of the classical case, is a CPTP map involving one 
neighborhood only: 



(14) 



where all the Vk(t) G ii(H m ) are quasi-local with respect to 
the neighborhood Afj. One of the reasons for focusing on this 
class of evolutions stems directly from applications: methods 
for implementing unitary evolutions, as well as related unital 
channels with the aid of some ancillary systems, are available 
in a number of diverse experimental settings. On the other 
hand, constructing arbitrary quantum channels is a more 
challenging task 091 . and can be generally done with good 
approximation only in the limit of fast control and/or short 
time scales [40|. The building block ( fi4] > can lead to different 
evolutions for the whole system, depending on neighborhood 
selection: 

• Random single interaction: at each time t one neighbor- 
hood Afj(t) is selected at random, j(t) being a single- 
valued random variable onto the neighborhood index set. 

• Cyclic single interaction: at each time t one neighborhood 
AfjU) is selected deterministically, for example periodi- 
cally cycling between the available j. 

• Random or cyclic asynchronous interactions: similar to 
the previous options, but a subset of several neighbor- 
hoods is selected at each time t. We can request the 
selected neighborhoods to be disjoint or not. This choice 
has an effect for implementation, but not for convergence 
properties of our algorithm, so we will not consider it 
further. 

• Synchronous interaction: all the available interactions are 
activated at each time, weighted by some qj > with 
Y^j qj = 1 in order to maintain a valid trace-preserving 
map: 

£( P ) = J2ij^Ap)- (15) 
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• Expected evolution: we study the evolution in expectation 
of the random interaction protocol which selects neigh- 
borhood Afj with probability qj at each t. Remarkably, 
the evolution to p t +i given p t then follows the same 
law ( fl5] > as the synchronous case. Note that convergence 
of the expected evolution to consensus does not guar- 
antee (at all) that a(ny) single evolution, determined by 
a realization of the random process {j(t)} t >o, would 
converge to consensus. Nevertheless, the statistics of any 
measurements performed at any time on the system will 
be exactly the same for (JT3J as for the associated random 
evolution. In this sense, convergence in expectation is 
indistinguishable from trajectory-wise convergence. 
The last two cases involve a time-independent map. Another 
time-independent map is obtained if we consider cyclic inter- 
actions of period T and we focus on the state at the end of 
every cycle: 



Pt+r = £c(ftt) = £m t ° ■ • ■ ° £jvi (fit)- 



(16) 



The consensus goal can now be specified formally. 

Let d(p a ,C) = inf peC \\p a ~ p\\, where C C £>(%) and 
|| • || is any p-norm on 05 (H). Given a sequence of channels 

{£t (•)}£()> define £t(po) = Pt = £t ° £t-i ° •• • £i(po)- 

Definition 6 (Asymptotic Consensus): A sequence of chan- 
nels {£ t (-)}^ , is said to asymptotically achieve aEC if 



lim d(£t(po),Ca-Bc) = 0, 



(17) 



for all initial states pq. 

The same definition holds for RSC, SSC, and erSMC by 



substituting the corresponding state sets in (17 1 



Definition 7 (Asymptotic Average Consensus): We say that 
the sequence of channels {£t(-)}t^o asymptotically achieves 
S -average aEC for some S G Sj('H m ) if it asymptotically 
achieves aEC and for all po, it holds: 



lim Tr(ap e (t)) 

t— ¥ OO 



= lim Tr(o-W 

t— J-OO 

= Tr(S> ) 



p(t)) = lim Tr(Sp(t)) 

t— >oo 

(18) 



for all I € {1, ... , m}. The same definition holds for ctSMC. 
We say that the sequence of channels {£*(•) It^o asymptoti- 
cally achieves S-average RSC ( resp. SSC) if it asymptotically 
achieves RSC (resp. SSC) and for S G S)(U m ) there exists a 
a G fi(H) such that (|T8j holds for all p . 

By expressing the action of quantum channels in the dual 
(Heisenberg) picture, it is possible to obtain a clear character- 
ization of the dynamics that satisfy ( fT8| . 

Proposition 6: Consider a sequence of CPTP channels 
{£t(-)}fZ Q , and call £ t = £t ° £t-i ° • • • ° £%■ The associated 
dynamics satisfies ([18) if and only if 



for 



S = lim £J{S) and lim £j(a 

t— >oo t— >oo 

1,2, ... ,m, where £[ — £\ o £} ~ 



S 



(19) 



of 



Proof: The conditions ( fT9] > clearly imply (TTHJ. On the 
other hand, if ( fT8] l holds for all p a , it is easy to obtain ( fT9| ) 
by duality, taking the limit inside the trace functional. ■ 
The first of the equalities in ( fl9] > holds in particular for the 
natural situation where £}(S) = S for all t. 



IV. A Gossip Algorithm for Quantum Consensus 

We now propose actual interactions that drive the quantum 
network to average consensus. As a building block, we focus 
on the interaction between two subsystems while the others 
remain unchanged; all neighborhood-activation options build 
on this elementary case, as explained above. 

A. Another viewpoint on the Classical Gossip Algorithm 

The standard linear consensus algorithm corresponds to 
(JS),© with f(x,y) = a(x — y). Its form with a single in- 
teraction activated at any time — also called gossip algorithm 
— is usually described as follows BP . At each iteration, a 
single edge (j, k) is selected from the set E(t) of available 
edges at that time. The associated agents move towards each 
other / their mean value, according to: 

Xj(t + 1) = Xj(t) + a(x k (t) - Xj(t)) 

= (l-f3)x 3 {t) + p x > (t)+ 2 Xk(t) 

x k (t + l) = x k (t) + a(Xj(t) - x k (t)) 

= {l-[3)x k (t) + l3 Xi{t)+ 2 Xh{t) 

x e (t + 1) = x t (t) for all I $ {j, k} , (20) 

where a G (0, 1) to have meaningful results^] and j3 = 2a. 
An alternative viewpoint on this behavior is that the inter- 
acting agents take a weighted average between two discrete 
operations: [keep your state] and [swap your state]; namely 

( Xj {t + 1), x k (t + 1)) = (1 - a) (Xj(t), x k (t)) 

+ a (x k (t), Xj (t)) (21) 
x e (t + 1) = x e (t) for all I £ {j, k} . 

This latter viewpoint turns out to have a natural quantum 
counterpart. Working with neighborhoods, one could also 
apply multi-agent permutations, e.g. 

(xj, x k , xi) t+1 =(1 - a- f3) (xj, x k , xi) t + a (x k , x h Xj) t 

+ P (xi, Xj, x k ) t . 

The following result (see e.g. BTI ') characterizes convergence 
to consensus with the gossip algorithm. While it is a known 
result, we nonetheless provide a proof that will be useful to 
our aim, i.e. proving the convergence of the quantum gossip 
algorithm, and makes our presentation more self-contained. 

Proposition 7: Consider G(V, E) an undirected graph that 
is connected, i.e. for any pair of vertices a, b G V, there exists 
a sequence of vertices vq = a, V\, V2, Vn-i, v n = b such 
that (v k _i,v k ) G E for all k — 0,1,..., n. If one step of 
the classical gossip algorithm ( |20] > is applied at each time, 
selecting the updated edge by cyclically running through all 
the edges of G(V, E), then the system exponentially converges 
to average consensus. Moreover, if the updated edge (j, k) is 
selected randomly according to a fixed probability distribution 
{qj. k }, with all qj k > 0, then asymptotic average consensus is 
ensured with probability one, in the sense that: for any 6,e > 
0, there exists a time T > such that 

P [\\x k (T) - x)\\ 2 > e\\x k (Q) - x\\ 2 ] <S, 

3 That is, the new states are an interpolation and not extrapolation between 
Xj(t) and Xk(t). 
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where P denotes the probability measure induced by the 
randomization, ||ir|| 2 = J2k x k x k and x = ^ l J2k x k( t ) = 
J2k x k{ty f° r an Y choice of the edges. 

Proof: We denote x T x = \\x\\ 2 for short and #E the 
number of edges in G(V,E). At any step of the gossip 
algorithm, W := ^ J2k,j IK ~ x jf = Efc \\ x k ~ x\\ 2 can 
only remain unchanged (if the two nodes of the selected edge 
have the same value) or decrease (as soon as an edge with 
different node values is selected). A Lyapunov argument on 
W then shows that the system must asymptotically converge 
to average consensus when the edges of a connected graph 
are selected in a cyclic way. Since the map associated to 
one full cycle of edge selections is linear and time-invariant, 
this convergence is exponential. For such convergence to be 
possible, there must exist some A > and integer M > such 
that W(T) < W(0)X if the edge choice between t = and 
t = T = M#E corresponds to M cycles of gossip iteration. 
When edges are selected randomly, any particular sequence 
of b consecutive edge selections has a probability greater than 
q b > to appear at least once during any time interval of 
length at least b, where q = min^ Mgg qj t k- In particular, if 
we target W(T) < eW(0) = X r W(0), we can say that there is 
a probability at least q rM # E to select r times a succession of 
M cyclic interactions between t = to and t = to + rM i^E. If 
this happens once, any preceding or following edge choice can 
only improve W (because of our first statement in this proof). 
We conclude by noting that over a time interval brM#E, 
there is then a probability < (1 — grM#Ey to have never 
selected r times a succession of M cyclic interactions, and 
thus potentially miss W(T) < eW(0); the probability that 
this happens can be made arbitrarily small by taking b (thus 
T) sufficiently large. □ 

B. Quantum Gossip Interactions 

Let us introduce a way to implement gossip-type interac- 
tions in a fully quantum way. In a controlled quantum net- 
work, one can typically engineer unitary transformations that 
implement the "identity" evolution and the swapping of two 
neighboring subsystem states; let us denote the latter operator 
by Ufj,k) f° r swapping subsystems j and k. By conditionally 
associating these two operations to orthogonal states |£f) and 
|£g) of an ancillary quantum two-level system and starting that 
auxiliary system in the state p 5 = (1 — a) +a|£s)(£sl> 
the joint state would evolve as: 

{l-a)p®\Zi)(Zi\+aU m pUl jtk) ®\Z s )(Zs\. 

The evolution of the composite system is thus described by 
a single unitary transformation^ Taking the partial trace 
over the ancillary system, we obtain as evolution for the 
quantum network a quantum gossip interaction implementing 
the quantum channel: 

p(t + 1) = SjMt)) = (l-a) p(t) + a U Utk) p(t)Ul k) , 

(22) 

4 The unitary transformation being given by I (g> 1 + Urj j.j Cg 1 

lfa)«sl- 



with a 6 (0,1). The optimal quasi-local mixing is obtained 
with a = 1/2. 

To develop our analysis, it will be convenient to introduce 
the graph Q associated to the multipartite system: its nodes 
1, ...,m correspond to the "physical" subsystems, the edge 
(j, k) is included if the subsystems j and k have a non-zero 
probability to interact. 

C. Convergence to Consensus 

We study convergence under three types of gossip dynam- 
ics: cyclic interaction, expectation of random interaction, and 
trajectory-wise for the random interaction. In all these cases, 
quantum gossip can be described by unital CPTP maps. We 
begin by recalling a characterization of the fixed points of such 
maps (see e.g. [19]). 

Proposition 8: Let {^}|Li the Rraus decomposition of a 
unital CP map £ (•) and define: 

A £ = {X e <8CH m ) \ [X,Vi] = 0Vi=l,...,K}. (23) 

Then X € Q3('H" 1 ) is a fixed point of £, i.e. £{X) = X, if 
and only if X £ Ag. □ 

This helps determine the set of fixed points for the CP maps 
of interest in quantum gossip. 

Lemma 1: Let Ufj^ denote the pairwise swap operation 
of subsystems (j, k) on H m . If the graph Q associated to the 
system is connected, then the set of fixed points of any CP 
unital map of the form 

£(X)= qo X+ £ q] . k U [m XUl jk) , 

U,k)eE (24) 

with qo + ^2 q i> k = 1 ' q ° ' ^' k ^ > 

coincides with the set of permutation-invariant operators. 

Proof: According to Proposition [8] above, the fixed points 
are the X satisfying XUfj^ = U(j^)X, or equivalently 
llh k ^XU(j.k) = X. The latter expresses that X is invariant 
with respect to pairwise swaps on all the graph edges. It is 
well known that sequences of pairwise swaps on the edges of 
a connected graph generate the full set of permutations on the 
set of nodes, and so we get the conclusion. □ 

The following lemma shows how the contribution of the 
identity, i.e. the trivial permutation, in the CP map plays a 
crucial role in the proof of convergence. 

Lemma 2: If a CP map £ admits an OSR with a term V\ — 
y/al with a > 0, then viewing it as a linear map on 23 ("H m ) 
its only modulus-one eigenvalue can be one. 

Proof: If £ is a CPTP map it is a contraction in trace 
norm lfl4ll . 1421 . so its eigenvalues belong to the closed 
unit disk. By virtue of the Kraus-Stinespring representation 
theorem (see e.g. l37l ). also T — jz^{£ — al) is CPTP and 
thus has eigenvalues fi k in the closed unit disk. Therefore the 
eigenvalues Afe = (1 — a)/ife + a of £ = (1 — a)F + al in 
fact belong to the circle of radius (1 — a) centered at a, which 
is strictly inside the unit circle except for a tangency point at 
1 e C, see Fig. [T] □ 
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Fig. 1, Delimitation of the (closed) domains for the eigenvalues {fJ-k} of T 
(blue) and {A^} of £ (red) [color online]. 

In other words, Lemma [2] excludes eigenvalues of unit norm 
different from +1, which are the ones that would cause limit 
cycles. 

By combining the above properties, we get the following 
convergence result for quantum gossip. It shows that S- 
average SSC can be attained for global operators that are the 
permutation-invariant average of local ones; this is similar to 
classical gossip, where distributed computation of the average 
of individual states actually gives access to the value of any 
linear permutation-invariant function of these states. 

Theorem 2: If the graph associated to possible interactions 
is connected, then the quantum gossip algorithm |22} ensures 
global convergence towards SSC: 

- deterministically, when the edges on which a gossip interac- 
tion occurs at a given time are selected by periodically cycling, 
in any predefined way, through the set of edges; 

- in expectation, when the edges on which a gossip interaction 
occurs at a given time are selected randomly from a fixed 
probability distribution {q j>k > 0| J2(j,k)eE = 1}: 

- with probability one on any trajectory, with the same edge- 
selection strategy of the previous point. Explicitly, there exists 
a state p* £ Cssc f° r which for any 5, e > 0, there exists a 
time T > such that 

P[Tr((p(T) - p*) 2 ) >e] <5. 

In any of the cases above, the state towards which the system 
converges is: 

P* = V U TPo Ul ■ (25) 
ml 

Furthermore, S'-average SSC is attained if and only if S G 
$)(H® m ) can be written, for some a € $j(H), in the form: 

S=-f> W (26) 

i 

Proof: First notice that all the operators in the OSR of 
the map |24) are self-adjoint. This implies that permutation- 
invariant observables S are fixed points for the associated dual 
map, and hence for the gossip iteration associated to any edge 
(j,k) and Vp : 

U U,k)SUl lk) = S =S> Tt[£ Jtk (p)S\ = Tr[p£j k (S)] = Tr[pS] 

(27) 



For the cyclic evolution map £c, we notice that all the 
simple two-subsystem swaps are still present with a weight 
different from zero in the OSR of the cyclic map ( fTo*| ), thanks 
to the presence of the identity in the OSR of each gossip 
interaction step. Therefore by Lemma[T]the fixed points are the 
permutation-invariant operators. Now consider the dynamics 
associated to £c as a linear, time-invariant map acting on 
the subspace of hermitian matrices. From Lemma [2] and the 
fact that the time-invariant linear map leaves D{H m ) invariant 
(excluding unstable Jordan blocks), we have that all the modes 
of the LTI system are asymptotically stable except those 
corresponding to the fixed-point set, namely the permutation- 
invariant set: every initial state converges to a fixed point Poo 
in this set. Thus the SSC set is globally asymptotically stable, 
and in fact exponentially stable since the map is linear. Let 
us now prove that Poo has the form |25|. For all permutation 
invariant X, from ( p7j ) we have that: 

Tr[X£ c (po)} = Tr[X Po ] V t. (28) 

Combining the latter with the fact that p,^ is permutation- 
invariant, that the set of all permutations is self-adjoint, and 
using (|3), we get for arbitrary Q g fj("H m ): 

Tr[Q Poo ] = Tr[Q—J2u„ Poo Ul] 
ml * — ' 

ml 

ml e - — ' 

7re<p 

ml L — ' 

jre<p 

This implies that indeed p x — p* as described in the 
statement. 

For the expectation of random evolution, the CPTP map £ 
is exactly of the form of Lemma [T] and the same reasoning 
can be repeated. 

For the random trajectory evolution, we repeat a proof 
similar to that of Proposition [7] Since £ for a single interaction 
is linear, self-adjoint, with eigenvalues in the closed unit 
disk, it is a contraction for the Frobenius norm distance 
Tr((pyi — ps) 2 ) between any two states pa,Pb € 2)("H m ). 
Indeed, £ has non-increasing orthonormal modes, so by writ- 
ing any operator X E Sj(H m ) in the modal basis we directly 
get Tr(£(X)t £(X)) < Tr(X^X); taking X = p A -p B yields 
the contraction. This is exactly analogous to the non-increasing 
Euclidean norm x T x = \\x\\ 2 under a classical consensus 
interaction with an undirected graph, and the related contrac- 
tion of \\xa — xb II 2 ■ Now taking in particular pa = p and 
Pb = p*, we get that the Frobenius distance from p to p* can 
never increase. Moreover, by transitivity of the permutation 
operators, i J2^ev U *P U l = h U nPo Ul = p* for 

any p along the trajectory of the gossip algorithm. Now given 
the convergence under cyclic evolution, there must exist some 
A < 1 and integer M > such that 

Tr((^(p)-p,) 2 )<ATr((p-p,) 2 ) 
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for any p for which ^ X^e<p U^pUl = p*. The proof 
then concludes along the same lines as Proposition [7] namely 
the probability to obtain an edge sequence which includes 
successions of M cyclic evolutions a sufficiently large number 
of times to have e-convergence, gets arbitrarily close to 1 if 
we wait long enough. 

Finally let us prove that we attain 5-average consensus 



if and only if S can be decomposed as in (26). We know 



from the first part of the proof that all permutation-invariant 
observables S are fixed points for the associated dual map £J. 
Then according to Proposition[6]we have .S-average consensus 
if and only if there exists a local observable a such that: 



lim £\{a^) = S 



for I = 1,2, ... ,m. Because of 
have that for every local operator 



and 



(29) 

251, by duality we 



lim E t (a 



1 

to! 



1 



5>». 

i=i 



(30) 



Therefore, we can achieve S'-average consensus if and only if 
S is a permutation invariant operator that can be decomposed 
as in d26|. □ 



Remark: This shows that the mean value of a (global) 
observable S = rprX)/ 17 > w i m arbitrary a, can be asymp- 
totically retrieved from the state of any single subsystem after 
having applied one of the quantum gossip algorithms. 

On the other hand, unlike for classical consensus, there are 
permutation-invariant operators that do not attain 5-average 
consensus, because they cannot be written in the form ( f2"6] >. 
This is the case among others if S is orthogonal to the linear 
span of all the local observables. For instance if S = af , 
given the orthogonal basis {<Jk}k=o,x,y,z for Q3("H), we have: 



Vfe G {0,x,y, z} and VI G {1, . . . ,m}. 



Therefore S cannot be written in the form Em, hence although 
S is conserved by the gossip algorithm, the latter cannot lead 
to 5-average consensus in the sense of Definition [7] 

As already mentioned, the convergence speed for the ran- 
dom case can be quite low, and a faster map would be obtained 
by effectively taking a mixture of all possible updates at each 
time. This can be attained by suitably selecting edges through 
an auxiliary quantum system, acting as an independent quan- 
tum "coin" (in the jargon of quantum random walks) at each 
time. This would represent a fully quantum implementation of 
a synchronous consensus update. 

D. Classical equivalent to observable consensus dynamics 

We next show how the quantum gossip algorithm ( p2) in 
fact implements in a quantum fashion the classical gossip as 
we restrict to ctEC. According to Definition [Tj a quantum state 
p belongs to C ct ec if: 

Tr[aWp] = ... =Tr[o-( m V]. (31) 

In view of this, it seems reasonable to attempt a convergence 
study of the algorithm (|22]i directly in terms of the evolution 



of the expectation values of the crW operators. This is not 
possible for arbitrary quantum evolutions, since a quantum 
state is far from fully specified by a single set of commuting 
observable expectations, and different states with the same 
expectation may lead to very different evolutions. However, 
our quantum gossip algorithm remarkably allows us to write 
a model for the average dynamics of the ctW in closed form. 
More precisely, let us define z#(t) :— Tr[£ t (p )a^} = 
Tr[pt<jW] . Note that for one subsystem swap Uu ^), we have: 



Tr[aWUu k) pU, 



(j,k)P u (j, k )i 



z t iftt{j,k} 
Zk ift = j 
z, if £ = k . 



(32) 



According to ( |32[ i and ( |22| ), the random gossip algorithm 
update yields, with probability qj,k, i.e. when the edge (j, k) 
is selected: 

{zj(t + l),z k {t + 1))= (1 - a)(zj(t),z k {t)) + a(z k (t),Zj(t)) 
z e {t+l)= Zi (t) for all I $ {j, k} . 

This last expression has exactly the same form as the clas- 
sical gossip algorithm ( |2"T| i. Therefore, Proposition [7] readily 
implies: 

Corollary 1: Under all the various edge selection strate- 
gies for quantum consensus algorithm (|22j>, the ze(t), £ — 
1, 2, m asymptotically converge towards the unique config- 
uration: 



^ m 

lim z t (t) = — E z fe(°) for a11 1 e {1,2,..., to}. 

m fc=i 

□ 

We remark that this only proves average a-Expectation 
Consensus of the quantum gossip algorithm, while our pre- 
vious Theorem [2] shows that the algorithm in fact ensures the 
stronger average Symmetric State Consensus. 



E. Gossip algorithm example 

In this section we briefly discuss the evolution induced 
by random quantum gossip interactions ( (221 on a four-qubit 
network whose associated graph is a patrj^J We observe its 
convergence toward average aEC, average RSC and average 
SSC. In particular we consider as a "target" global observable: 



S = 



i(o?> + (33) 



Let the initial state be: 



p = |1,0,1,0)(1,0,1,0|, 



(34) 



which is pure, and does not satisfy any of the consensus 
definitions provided in Section [TT] 

5 I.e. the available neighborhoods, labeling the subsystems as {1,2,3,4}, 
are {1,2}, {2,3} and {3,4}. 
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By Theorem [2] we have that the state asymptotically converges 
to: 

1 



Poo = lim p(t) 

t — yoo 



3! ^ 



(11,1,0,0X1,1,0,01 

h|l,0,0,l)(l,0,0,l| 
h|0,l,0,l)(0,l,0,l| 



|1,0,1,0}(1,0,1,0| 

|0,1,1,0}<0,1,1,0| 
|0,0,1,1}<0,0,1,1|). 



(35) 



This expression is clearly invariant under all the subsystem 
permutations, i.e. is in SSC, and therefore also in RSC 
and ctEC for all a. The expectation value of S is preserved at 
any step, and by Theorem [2] the algorithm drives the system 
to S'-average consensus, with a = a z . 

However, p^ is not in crSMC for any a ^ al. Indeed, 
according to Proposition |4] p,^ is in crSMC if and only if 
Trlpoollsym] = 1. Now let {IIi}f =1 denote the orthonormal 
rank-one projectors in ( |35| ) and define the orthonormal projec 
tor fl = £\ IT, such that Poo = § S 
get 



1 n, = in . We then 



Tr [poo n s y m ] 



1 6 



= ^Tr[n sym n]. 



(36) 



This last expression is equal to 1 if and only if Tr[n sym il] = 
6. However, excluding the trivial case a = al, for qubit 
networks H sym is always a two dimensional projector, so 
Tr[n sym II] < 2. Hence p^ cannot be in crSMC for any non- 
trivial a. 

Figure [2] shows the evolution of the expectation values 
of the local and of the global observables related to a z as 
the iterations proceed for one run. The edges are selected at 
random with uniform probability, and the mixing parameter 
a is taken to be 1/2. With this particular choice, the reduced 
density operators of two subsystems that have just interacted 
are equal; this explains why a maximum of three points 
are visible on the graph at any time. The plot shows that 
asymptotically the expectation of the local observables a z tend 
to the expectation value of the global observable S, while the 
expectation value of S is preserved at each step. 



1.0-QOOO 



0.5 - 



OOOOOOOO 



OOO 



♦ Tr[o- z 'p;] 
□ Tr[cr ; 2 p/] 
Tr[cr ; 3 p,] 

O Tr[cr- 4 p,] 
Tr[Sp,] 



0.0 - 



n 

««♦♦♦' 



-0.5 - nan 



10 



20 



30 



40 



Fig. 2. Evolution toward cr-Expectation Consensus for a four-qubit network 
arranged in a path graph. 



(a quantum network). We propose and analyze a quantum 
gossip-type algorithm that asymptotically prepares symmetric- 
state consensus states while preserving the expectation of any 
permutation invariant observable. 

A number of questions remain open. Among these, we 
believe that it would be particularly interesting to further ex- 
plore the link between single <7-measurement consensus states 
and entangled states, and to determine if, and under which 
conditions, it is possible to achieve this type of consensus 
with a distributed algorithm. This could potentially lead to a 
class of algorithms that prepare entangled states in a robust 
and distributed way. Another interesting point is to assess 
the potential of devising continuous-time quantum consensus 
algorithms. This could build on some sort of "continuous 
swapping" Hamiltonian dynamics and lead to connections with 
physically relevant many-body Hamiltonians and dynamics. 

Lastly, let us remark that in this paper we proposed a 
quantum algorithm in which the gossip-type interactions are 
selected in a classical way. The potential advantage of a 
fully quantum implementation, along with its connection to 
quantum random walks and Markov chain mixing properties 
fPTl . |[l"8l . is definitely worth further investigation. 



V. Conclusions and Research Directions 

In this paper we develop a general framework for posing 
and studying consensus problems in the quantum domain. In 
particular we provide various operationally-motivated general- 
izations of a "consensus state" to quantum systems - namely 
(T-expectation consensus, reduced state consensus, symmetric 
state consensus, and single er-measurement consensus - and 
establish their hierarchy. We highlight at each step the sym- 
metry considerations underlying the results, making explicit 
connection with the usual multi-agent consensus problem. 
The developments could be adapted in particular to obtain 
a "consensus on probabilities" framework for classical sys- 
tems. With respect to the existing work on non-commutative 
consensus (34], our approach follows the analogy with the 
classical setting as closely as possible, maintaining an op- 
erational viewpoint and working with a multipartite system 



Appendix 

A. Description of quantum systems and notations 

1) Quantum systems basics: This paper considers finite- 
dimensional quantum systems. Their mathematical description 
starts by considering a finite dimensional complex Hilbert 
space T~L ~ C d . The (Dirac's) notation \ip) denotes an element 
of T~L (called a kef), while = is used for its dual (a 
bra), and (ip\ip) for the associated inner product. We denote 
the set of linear operators on "H by 23 ("H). The adjoint operator 
E Q3("H) of an operator X 6 25('H) is the unique operator 
that satisfies (X|^)) t \x) = (V>| (* f lx» for all |</>>, \x) G U. 
We then denote S)(H) the subset of 25(H) of self-adjoint 
operators, and il('H) C 23(H) the subset of unitary operators. 
The natural inner product in 23(7^) is the Hilbert-Schmidt 
product (X, Y) — Tr(A^F), where Tr is the usual trace 
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functional (which is canonically defined in a finite dimensional 
setting). We denote by / the identity operator. Working in 
a finite dimensional setting, we often consider vectors and 
operators as represented by complex matrices of suitable 
dimensions: £ H ~ C d are represented by column vectors, 
so ((f>\ £ "H f ~ C d are row vectors; X G 93(H) ~ C dxd 
are d x d complex matrices, the adjoint X> is the transpose 
conjugate of X, self-adjoint and unitary properties carry over 
to the associated matrices. 

In statistical quantum theory, the state of a quantum system 
is represented by a density operator p, that is any self-adjoint, 
positive semi-definite operator with trace one. We denote the 
convex set of these operators (the state space) by D(H). The 
extreme points of this set, namely the rank-one operators p = 
with € H and = 1, are called pure states. 

A projective (or von Neumann) observation, or measure- 
ment, of a quantum system is characterized by a so-called 
observable, that is a self-adjoint operator a £ $j(H), see 
e.g. [36). Its spectral decomposition u = 2~2j s j^j determines 
the possible outcomes {sj} of the measurement, and the 
projectors Hj determine the associated state of the quan- 
tum system after the stochastic measurement: having state 
p before the measurement, the latter' s outcome will be Sj 
with probability P p (IIj) = Tr(ILjp) =: p 3 ; and if the j-th 
outcome is measured, then the state after the measurement 
is p\j = HjpHj / pj. The probability to observe s' k in a 
subsequent measurement of a' = s' k H' k , where the H' k 
do not necessarily commute with the Tlj, is then: 

p p |.(n' fc ) = Tr(n / fc n jP n j )/p i . 

From this it follows that the probability of observing the 
ordered sequence of two events first Sj, then s' k , given the 
initial p, is 

F p (u j ,u' k ) = Tr(n' k n jP n j ). 

If Hj and H' k do not commute, a different ordering in a se- 
quence of measurements can change the resulting probability. 
If IL, and H' k do commute, and only then, the joint probability 
of observing Sj , s k is independent of the measurement order 
for all p, and simplifies to 

p p (n„n' fe ) = Tr(n' fe n jP ). 

2) Multipartite systems and partial trace: For simplicity, 
we present the interaction of two quantum systems; the case 
of n > 2 systems is easily obtained by iteration. If two 
quantum systems, with associated Hilbert spaces Hi and H 2 
respectively, are taken together to form a larger bipartite 
quantum system, the Hilbert space Hi : 2 associated to the com- 
posite quantum system is the tensor product of the individual 
quantum subsystem Hilbert spaces, Hi <S> H?.- 

Let {|V'fe)}fcLi an d {|^)}f=i b e orthonormal bases for H\ 
and H2 respectively, then an orthonormal basis for Hi y 2 can 
be written as: 

im ® , (37) 

from which we get that dim(H 12 ) = dim (Hi) dim (H2) = 
d\d2- We use the short notation \i/},4>) := \ip) ® \<fr) for any 
\ip) € Hi and \<fr) £ H2- The composite Hilbert space is 



naturally endowed with the inner-product (1*1, U2\vi, 1)2) ■= 
(ui\vi) (u2\v2) ■ A representation and basis for operators in 
2K^i,2) is derived from its vector counterpart in the standard 
way. In particular, given two operators Xi £ 23 and X 2 G 
Q5(?^2), one can define Xi ® X2 £ %$(Hi,2) as the linear 
operator such that V|ui) € Hi , \u2) £ H2- 

Xi®X 2 {\ui)®\u 2 ))= Xi\ui)®X 2 \u 2 ) • (38) 

If two operators are in the form Xi an d I\<g)X2, i.e. they 
act non-trivially only on different parts of the multipartite sys- 
tem, then they commute for any X\ and X2- It is worth noting 
that in matrix representation, the tensor product corresponds 
to the Kronecker product. 

The partial trace over Hi is a linear map: 

Tr Hl : B(Hi ® H 2 ) — > B{H 2 ), (39) 

such that, for any X12 £ B(Hi,2) and any X2 £ it 
holds that: 

Tr[Tr^ [X 1>2 ]X 2 ] = Tv[X h2 (2i ® X 2 )]. (40) 

If now {IV'fellfcLi an< i {|^)}f=i are orthonormal bases for Hi 
and H 2 respectively, the partial trace over Hi can be written 
as: 

Tr Hl lXi, 2 ]=J2(^®&\Xi,2\^k®&M)(&\. (41) 

k,l,i 

The partial trace over H2 writes in a similar fashion. 

The tensor product structure together with the superposition 
principle enrich quantum theory with the phenomenon of 
entanglement. Namely, if some vectors in |£) £ H1.2 - like 
the above-defined basis vectors - can be factorized as 

10 = m ® 10) (42) 

for some \ip) £ Hi and \<p) £ H2, there exist many vectors 
in H1.2 which cannot be written like ( |42| >; these are called 
entangled. 

In the density operator formalism, a state p £ ?D(H± t 2) is 
called entangled if it cannot be written as a convex combina- 
tion of factorized operators, i.e.: 

P ^ Pl Pi®Pi Wkh Pi G aIld Pi G ®C^2) ■ 

(43) 

Intuitively, this means that an entangled state contains some 
specifically quantum correlation so that it cannot be separated 
into subsystem states (conditioned by classical probability 
correlations). The only individual characterization of a sub- 
system that we can give is its "expected" density operator 
if information about all other subsystems is ignored. This 
reduced density operator e.g. for subsystem 1 is obtained by 
taking the partial trace over H 2 of the overall state p: 

Pi = Tr« 2 [p] . 
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B. Distinguishing SSC from RSC by measurements 

Accumulating statistics about measurement outcomes at 
each subsystem separately, allows in principle to detect a pure 
state for which SSC and RSC are equivalent (see Proposition 
BJ. Indeed, the knowledge of all local measurement statistics 
is equivalent to knowing the reduced state. 

Proposition 9: Except for the case of Proposition [3] SSC 
can only be distinguished from RSC by inspecting correlations 
between measurement outcomes at different subsystems. 

Proof: The statement builds on the standard fact that the 
statistics of a local observable o\ %a%® ...®a m only depend on 
reduced states pi, p 2 , Pm- So repeated local measurements 
can, at their best, fully characterize the p k . Checking RSC, 
i.e. that these p k are all equal, is thus straightforward. On the 
other hand, reduced states p k are the best that can be extracted 
by local measurements in trying to distinguish RSC from SSC 
states. If pi = p2 = ...=: p have rank one, we have the special 
case that is always SSC. If instead p has rank at least 2, we 
can write it as p = p\R\ + P2R2 where Ri, R2 <E 'B('H), 
Pi, P2 are positive scalars, R2 is positive semidefinite, and 
Ri is a projector on a 2-dimensional subspace V 2 . Consider 
Ri = |e 1 }(e 1 | + |e 2 }<e 2 | = |/i}</i| + |/ 2 }</ 2 |, where |e x >, |e 2 ) 
and \fi), I/2) are two orthonormal bases for V 2 with (ei|/i) ^ 
{0, 1}. Now the reconstructed p could equally well reflect the 
state 

P = P® m , 

which is SSC, or e.g. a state of the form: 

p = P2 Rf m 

+ MleiM) + |e 2 )|/2))(|e 1 )|/i) + |e 2 )|/ 2 ))t ® R® (m " 2) , 

where the first two subsystems are entangled. This state is not 
SSC, even for m = 2. Thus the local knowledge of p does not 
allow to distinguish if the state is SSC or not. □ 



n s y m II S y m II S y m p 



p, 



where lT, ym is associated to a' . Since H := 



nsymn s y m ii s y m 



and p both are self-adjoint positive semidefinite, the only way 
to have Hp = p 7^ is if H has at least one eigenvalue 
> 1. Now take in particular a — ^ fe k\x k ){x k \ and a' = 



(thus the 



E k k\ Pk ){p k \, with p k = ^ ^=0 e^ n \ X] ) 
|p/c)-basis is related to the |xfc)-basis by Fourier transform). 
A few computations show that H then has all eigenvalues 
< 1, except for n = 2 that is the case of two qbits. For the 
latter particular case, one can prove the property by showing 
e.g. that there is no state which would satisfy crSMC for all 
a € {cr x ,a y ,a z }. □ 
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C. Proof of Proposition [7] 

The definition of trSMC involves T^II^II^p). which [17 
takes the partial trace over the state of all subsystems except 
the pair {k,£}. So we can effectively discard all but two 
subsystems, and show without loss of generality that it is 
impossible to make crSMC hold for all a on two subsystems 
k = 1, £ = 2. In Proposition [4] we say that trSMC for a 
particular a requires U sym p = p with II sym = . IT,® m , and 
{ITj} the spectral projectors associated to a. So if ctSMC has 
to hold for both a and a' , we must have in particular 
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